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We rewrite the exact expression for the finite temperature two-point correlation function for the 
magnetization as a partition function of some field theory. This removes singularities and provides 
a convenient form to develop a virial expansion (the expansion in powers of soliton density) . 



To calculate correlation functions in strongly corre- 
lated systems is not an easy task, even if the corre- 
sponding models happen to be integrable. For models 
with dynamically generated spectral gaps the most pow- 
erful technique is the formfactor approach pioneered by 
Karowski et. al. (H and perfected by Smirnov This 
approach works wonderfully for zero temperature, but 
encounters difficulties at T ^ 0. These difficulties are 
related to singularities in the operator matrix elements 
(formfactors). These singularities exist for operators non- 
local with respect to solitons, they originate from forward 
scattering processes and their treatment requires careful 
infrared regularization. Despite long efforts a correct reg- 
ularization has not been yet found. However, for models 
of free fermions (such as the XY model or the Quan- 
tum Ising model), there are alternative means to calcu- 
late the correlation functions which allow to bypass the 
above problems. These alternative approaches include 
the determinant representation of the correlation func- 
tions 0|,f4] and the semiclassical method (which may 
have much wider apphcation, see For these re- 

sults to have a greater use one has to establish their re- 
lationship with the formfactor approach. A step in this 
direction was made in |3| where the semiclassical results 
[llil^iQ were reproduced by summing up the leading sin- 
gularities in the formfactor expansion. Such summation 
was restricted to the leading order in the soliton density 
n ^ exp(— A//T) {M is the spectral gap). In this paper 
we describe a formfactor-based representation for corre- 
lation functions which, though in its present form is valid 
only for models of free fermions, is rather suggestive and 
may give rise to useful generalizations in the future. For 
the Quantum Ising (QI) model, which is the main ob- 
ject of this paper, this procedure naturally gives rise to 
a virial expansion of the dynamical spin susceptibility. 

The QI model Hamiltonian is 

H = 5]hJ«+i-H/i<] (1) 



where s^,s^ are the Pauli matrices. By the Jordan- 
Wigner transformation this Hamiltonian can be trans- 



formed into the Hamiltonian of non-interacting fermions: 

H = Y, <P)F+Fp; e{p) = ^ {J - hf + AJ h s\n\p / 2) .{i) 
p 

Hamiltonian |^ possesses a property of self-duality: the 
transformation = ^^J' = «+i 

preserves both the commutation relations and the form 
of the Hamiltonian: 

n 

As follows from the form of the dispersion law h = J 
is the critical point. It is easy to see that at T = 
it separates the regions where (cr^) 7^ 0(J > h) and 
7^ ^("^ ^ Operators a'^ and /i^ are re- 

spectively called order and disorder parameter operators. 
The duality allows to study the correlation functions at 
one side of the transition only. For instance, a correla- 
tion function oi cr^ aX J > h coincides with the correlation 
function for /i^ with J, h interchanged. 

Though one cannot observe a single fcrmion, the 
fermionic statistics can be indeirectly tested by measur- 
ing the correlation functions of . The operator cr^ (as 
well as n^) is local in fermions: 

fe 

-f{k) = ^l + k/e{k). (4) 

Its finite temperature correlation functions are just polar- 
ization loops; they clearly contain the Fermi distribution 
functions of the fermion fields F. Since the order {a^) 
and disorder parameter (^^) fields are nonlocal in terms 
of fermions, their correlation functions are more compli- 
cated. 

The duality allows us to consider h > J 01 h < J phase 
only. In this paper we will study the two-point correla- 
tion functions of the and /i^ operators in the "ordered" 
phase h < J. For technical reasons it will be convenient 
to work in the limit | J — /i| << J when the spectral gap 
is much smaller than the bandwidth and one can formu- 
late a continuous description. In this limit the excitation 
spectrum is relativistic e{p) — ^ (?p^ + Af ^ , (? = Jh. 
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Energy and momentum of a quasi-particle are conve- 
niently parameterized by a rapidity, 0, {cp = M sinh 0) . 
Then the eigenstates of Hamihonian ^ are labeled by 
sets of rapidities, {9i}, such that the energy and momen- 
tum of the system are equal to 



E = M^cosh^ 



P 



c^^M^ sinh 6*,. (5) 



Below we set c — 1. 

A convenient finite temperature expression for the two 
point correlation functions of a and /i was derived by 
Bourgij and LisovyyQ. This expression for the Matsub- 
ara time correlation function is manifestly free of singu- 
larities and has the following form: 



(a(T,x)a(0,0)} = CMi/4e-l"l^(^) x 



E 



(2iV)! 



E n 



n 



(6) 
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where r is imaginary time. The same expression holds 
for fi^, but with 2N replac ed by 2N +1. q ^ 2TTTm 
(m integer), and e{q) = y^m^ + q^. The term in the 
exponent is {/3 = 1/T) 



■nil) = 



2e{q) 



Ax 



, „ . , „ lncoth[/3e(a;)/2], (7) 



and 



A(T) 



dp 



ln{coth[/?e(p)/2]} 



(8) 



The symmetry breaking transition at T = leads to a fi- 
nite magnetization, (a) = ±[CM^/^]^/^. This is reflected 
in the zeroth order term in Eq. 0. 




FIG. 1: A graphic representation of Eq.®. The ellipsises 
are formfactors of a operator. Lines with left arrows are 
/'•"'"■' (^) exp[ite(S) + ixp{9), lines with the right arrows are 
/(-'(0)exp[-it6(0)-ia;p(e)]. 

In % Eq. ® was rewritten in the form which allowed 
an analytic continuation for real time. The iV-th term in 
the square brackets of Eq. ^ is given as 0| 



27V 

E 



n\{2N 



i=l 



2N- 



n 



-\xp. 



n 



i>k 



tanh'^ [ 



-^/c)/2]n,>ptanh2[(0^-0;)/2] 



n,.ptanh^ 



0' +i0)/2] 



(9) 



where are numbers of particles and antiparticles. Now we rear- 

, , range the double sum 151) in such a way that we first sum 

»)' '(e) -77'~'(e) . . . 

f{+)(Q.] — _^ f^^Uf)-] — — ^ (^0) terms which contain a fixed difference between num- 

[e/3«(e) _!]'•' ^ [1 „ e-0<o)Y bers of particles and antiparticles 2N - 2n = 2k. Then 

(±)(n\ i-^sinh6' f°° da; ln{coth[/3e(x)/2]} such term in Eq.Q can be represented as an integral of 

^ ^ > ^ J ^ ^2 _ ]\/[2 sinh^(6' ± iO) correlation function of a Gaussian field theory: 

and p — M sinh^?, e — M cosh 6. In Eq.©, n and N — n 



/N-k ,„ N+k 1/1/ N-k N+k 

TT ^f(+)mAer^.+i-p. TT ^ fi-) mAe-^-i-n /T] e^'^ie.+ia) TT g-i*(e;-ia)g2fei$(oo)^n 
1=1 j=i % j 



where and a >> oq — s- 0. Ea. (|ll|l represents one term in the 

perturbative expansion of the partition function of the 

($(0i)$(02)) = Goieu) = - In [tanh2(^^l2) + 4] (12) 
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theory with the action 



d9 
27rao 



v[0,m] 



(13) 



Y _ j(+)j-^-jgTMcoshe+ia:Msinhegi*(e+ia) ^ 

^ j(-) j-^-jg-rM coshe-ixM sinh eg-i$(e-ia) {lA) 

In this theory x, t are external parameters. The field $ 
lives on an infinite line in 9 space. It is easy to see that 
the entire correlation function can be written as 



(15) 



fc— — oo 



and {■■■)v stands for averaging with action (|13|l . Already 
this form leads to significant simplifications. Namely, 
the partition function can be written as an exponent of 
the free energy Zq = exp[— F(a;, t)]. The latter one is 
represented as a sum of the cumulants: 



oo , r ^ Aft 



N 



- (iV!)^ 



N 



N 



where the subscript stands for averaging with V = 
and 



oo „N+k j„ N — k 1/1/ 

(( n gi^^^'+i'^) n g-i*(e^.-ia)g2ifc$(fl)^^^ (^7) 
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All transformations so far have been exact. Now we 
would like to concentrate on the causal Green's func- 
tions. To obtain them one has to replace r with it in 
©. Assuming that T <C M, we consider the region of 
frequencies \ui\ < 2M, where the only terms of the expan- 
sion @ contributing to the spectral function are those 



which contain equal number of particles and antiparti- 
cles. As we shall demonstrate, all formf actor singularities 
are contained in the first term of expansion (|16() . Since 
/(+)| exp(-/3Af) << 1, |/(~^| 1, this expansion is 
in powers of soliton density exp(— /3Af). The first term 
is given by 



47r^ J tanh (6'i2 + i0)/2 

-1^ I dM./(+)(^ + .)/(-)(^-.) ^"P["^(t""^^t"'°'^^^^ (18) 

9{t-\x\) \ - [ d9f'^+\9)f<--\9)M[{tsmh9 + xcosh9)\] + -g[9 = tanh-i(x/i)] \ 

I Jta.nhe<-\x\/t J 



where g{9) ^ f^+\9) - f'^-\9) and, as follows from ^ 
/(+)(6))/(-)(6/) ^ {4cosh^[/3e(6')/2]}-i.Here we assumed 
that Mt,{Ttf/'^ > [1 - (a;/i)2]-i/2. The higher order 
cumulants contain higher powers of exp(— /3M) and also 
do not contain positive powers of t. This justifies keeping 



the exact distribution function in the real part of H18(l . 
This cannot be done for the imaginary part, since the 
second cumulant gives a time independent contribution ~ 
exp(— 2/3i\f ). Within these limits we obtain the following 
results: 
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{a{x,t)a{0,0))T = CMi/''6'(t) exp(-(5A|x|)exp 



1 

47r 



dp 



\tv{p) - x\ 
cosh2[/3e(p)/2] 



4i 



{a{x, t)a{0, 0)}t = CM^/*e{t) exp(-A|a;|), |a;| > t 
{fiix,t)^iiO,0))T = {fi{x,t)fi{0,0))T=oHx,t)aiO,0))T 



eM-0M/^l-{x/t)^]j,t> \x\ 

(19) 
(20) 



SA 



dp 



lncoth(/3e/2) 



2cosh^(/3e/2) 



exp(-3/3M) 



The imaginary part of (|19|l in the time-Uke domain t > |x| 
reflects a quantum nature of the excitations. For T — 
such imaginary part was first found in jlll | . In the leading 
order in exp{—(3M) Ea. H2U|) coincides with the one found 
in 0. At x = we have 



(cr(x = 0,<> 0)o-(0,0)) 



(21) 



CM^/^ exp 



^t/ro + -e-P'' + O (e-2^*^) 

TT 



where Tq = ^ e^^^+i • ^^'^•^'3 single solitons are not di- 
rectly observable, it is rather interesting to note that tq 
contains the distribution function of a single soliton. It 
would be very interesting to see what happens in corre- 
lation functions in models containing particles with frac- 
tional statistics. 

Since our approach shares certain common features 
with the Fredholm determinant representation intro- 
duced by Korepin et. al. 0,11, [13,111, we feel obliged 
to comment on the subject. The main diff'erence is that 
we do not represent the correlation functions as deter- 
minants though in certain limits this is possible. For in- 
stance, if one adopts the nonrelativistic limit 6* << 1 in 
action (|13|l . it can be fermionized and rewritten as a the- 
ory of free fermions. Then by integrating over fermions 
one obtains the determinant representation. However, 
we would like to point out that such representation is 
not a goal in itself. By representing correlation functions 
as partition functions of some field theory one already 
achieves a lot since now one can concentrate on connected 
diagrams where it is easier to keep track of sing ularities. 
It is possible that acting along the lines of [lj| one can 
obtain such representations for interacting models. 

We also would like to warn against the direct compari- 
son of Ea. (|19|l with a similar equation for the ((T^cr+) cor- 
relation function in the XY model obtained in 3j. This 
warning is necessary because the XY model in magnetic 
field is rather similar to the QI model; the similarity in- 
creases when the magnetic field exceeds the band width 
so that the ground state becomes ferromagnetic. How- 
ever, as was explained to us by Korepin (private commu- 
nication), the formulae of '3| were obtained for the case of 
weak magnetic field the XY model spectrum is a gapless. 



which explains the difference in the final results. 
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